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Abstract - We review recent theoretical and experimental results demonstrating the interac- 
tion between force fields and density inhomogeneities as they arise in premixed turbulent 
flames. In such flames the density fluctuates between two levels, the high density in reactants 
p r and the low density in products p p with the ratio p r /p p on the order of five to ten in flows 
of applied interest. The force fields in such flames arise from the mean pressure drop across 
the flame or from the Reynolds shear stresses in tangential flames with constrained stream- 
lines. The consequence of the interaction is nongradient turbulent transport, countergradient 
in the direction normal to the flame and nongradient in the tangential direction. The theoreti- 
cal basis for these results, the presently available experimental support therefor and the impli- 
cations for other variable density turbulent flows are discussed. 


137 



INTRODUCTION 


The phenomenology of turbulent flows is largely based on gradient transport assumptions proposed over a 
century ago by de St. Venant and/or Boussinesq. Even in the recently exploited second moment methods 
of analysis, methods which involve large systems of partial differential equations calling for considerable 
computational effort, closure is achieved by employing gradient models to eliminate third moment and 
other quantities. The extension to turbulent flows with variable density of the various gradient models 
carefully developed and validated for constant density turbulence has unfortunately been casually under- 
taken with the consequence that much additional research is needed before the phenomenology is well 
founded for high speed turbulent boundary layers and jets and for turbulent reacting flows, both high and 
low speed. That new significant processes may be operative in turbulent flows with variable density is 
suggested by the recent findings in premixed turbulent flames. It is our purpose to review these findings 
and to discuss their implications. 

Exposition is facilitated if we consider a normal premixed turbulent flame as shown schematically 
in Fig. 1. Cold reactants, a metastable mixture of fuel and oxidizer, with a mean velocity u 0 , a density 
p r , an intensity of velocity fluctuations characterized by a turbulent kinetic energy k a and a length scale 
of the large eddies l Q are consumed within the flame and exit as hot products with these quantities 
changes to (1 + x) u a , p p , k^ and /„ respectively. Here x is a heat release parameter with values of practi- 
cal interest from five to ten. Figure 1 may be considered an idealization of the flames which occur in 
internal combustion engines, various propulsion systems involving prevaporized fuels and industrial 
accidents. 

The object of a theory for such flames is the prediction of the turbulent flame speed u 0 and of the 
flame structure. The development of such a theory involves two related considerations; one concerns the 
thermochemistry of the flow, the description of the state of the gas, its density, temperature and composi- 
tion, while the second relates to fluid mechanics. These two aspects are the aerothermochemistry of the 
flow. 
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Various parameters can be used to characterize these flames. One system is proposed by Abraham 
et al. (1985) and involves a Damkohler number Da A = (A Ju ') (S L /b L ) where A is a large eddy scale, i.e., 
l a introduced earlier, u' ~ (3/2 k a ) v2 is a characteristic intensity of the velocity fluctuations and S L and 
81 are the speed and thickness respectively of an unstrained laminar flame in the chemical system under 
consideration. A second parameter is a turbulent Reynolds number R A = u' A/v where v is a representa- 
tive kinematic viscosity.* A pair of values, Da A and R A , determines the ratio of two velocities u /Si and 
of two lengths L k / 8 l where L K is a characteristic Kolomogoroff length. For premixed turbulent flames 
within internal combustion engines Abraham et al. (1985) determine that 1 < ( L k / 8 l ) < 10 2 . The impli- 
cation of this finding is that chemical reaction at the molecular level takes place in thin surfaces, laminar 
flamelets, whose structure is dominated by laminar transport and whose motion is determined by the tur- 
bulent velocity field. Similar considerations applied to other premixed turbulent flames indicate that in 
many, indeed most, applications of applied interest this laminar flamelet description prevails. 


The Bray-Moss Model and Some Consequences 

Although not restricted to flows involving such flamelets, the Bray-Moss (BM) model simply and with an 
accuracy suitable for many purposes describes the thermochemistry of premixed turbulent flows (cf. Bray 

and Moss 1975). According to this description the instantaneous value of a progress variable c(x,t) 

which has the values zero in reactants, unity in products and intermediate values in gas possibly undergo- 
ing chemical reaction determines the entire state of the gas. Thus, for example. 


= X = 1+Tc (l) 

p r 1 + x c T r 

where x is most readily determined from the specialization of the second of these equations to products, 
i.e., from 



* It should be kept in mind that V can vary by one to two decades within a flame and thus that there is some ambiguity as to its appropriate value 
for determining R A - 
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In averaging the equations for variable density tuibulent flows it is useful to employ Favre or mass 
averaging; to illustrate the notion and notation involved consider the i-th velocity component, namely 


Ui(X,t) = 


p U; 

-=-(x) + U "(x,t ) = u(x) +u "(x,t) 

p - ~ - 


All variables except the density and pressure are averaged in this fashion. The great advantage of Favre 
averaging is that the conservation equations for variable density flows closely resemble those for constant 
density with exceptions which, as we shall see, indicate special effects associated with the variability of 
the density. Despite this advantage it should not be assumed without support from experiment that the 
closure methods applicable to constant density flows can be carried over to variable density cases without 
modification. Furthermore, although significant only when low turbulence Reynolds numbers must be 
taken into account, Favre averaging is not adaptable without approximation to the molecular terms in the 
conservation equations. The definitive discussion of Favre averaging applied to turbulent combustion is 
given by Jones (1980) while such averaging is emphasized in the monograph on turbulent reacting flows 
by Libby and Williams (1980). Rubesin and cowoikers at the Ames Research Center of NASA apply 
Favre averaging for the analysis of high speed turbulent boundary layers (cf. Wilcox and Rubesin 1980). 

When Favre averaging is applied to Eqs. (1), there results 


P = 1 

Pr 1 + T C 


T 

~=\+XC 

T r 


( 2 ) 


A further important feature of the BM model is an approximation for the probability density func- 
tion (pdf) of the progress variable. In general for a statistically stationary flow 

P (c \x) = aix) 5(c) + p(x) 5(1 -c) + ^x)f(c-,x) ( 3 ) 

Now if 

r 

J dc f(c ;x) = 1 

cr 

i.e., if the pdf describing the interior values of the progress variable is normalized, then 
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a(x) + p(x) + y(x) = 1 


(4) 



The essential approximation in the BM-model is that y« 1, an approximation which is satisfied when 
chemical reaction occurs in laminar flamelets as well as in other flow structures, e.g., reaction zones with 
5/, > L k . The implication of this approximation is that the temperature as measured by an idealized ther- 
mometer within a premixed turbulent flame possesses essentially two values, T r within reactants and T p 
within products, and that flamelet passages cause the switch from one level to the other.* 

An important consequence of the y « 1 approximation is exposed if c (x) is calculated from Eq. (3). 
There results 


a(x) = 


1 — c 

1 + TC 


P(*) = 


(1 +T)C 
1+TC 


(5) 


where Eq. (4) is used to calculate ct(x). We thus see that the strengths of the delta functions is simply 


related to the mean value of the progress variable which must be considered one of the principal depen- 
dent variables in a theory of premixed turbulent flames. 


For the present discussion it is useful to consider the intensity of the density fluctuations; from 
Eqs. (1), (3) and (5) we can compute 


P /2 ,2 c (1-g) 

T i + T 


( 6 ) 


We thus see that the maximum in the relative intensity of the density fluctuations occurs within the flame 
where c = 1/2 and has a value t 2 /4(1 +t) = x/4 for % » 1. Thus for degrees of heat release of applied 
interest we must expect intense density fluctuations. 

Although not essential for the present discussion, it is worth noting that the mean rate of creation of 
product vv(jc) in the BM model involves the product y w max where w max is the maximum rate of creation 


of product, a quantity determined by the chemical kinetics of the system under consideration. In the BM 
model this product involves a vanishingly small term multiplied by an indefinitely large term and is thus 

* Extensive studies of the statistics of two valued functions with specific reference to premixed turbulent flames have been carried out in order to 
develop a model for the mean rate of creation of product (cf. Bray et al. 1987). 
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indeterminant. As a consequence a separate model for w(x) is needed; progress in this direction is 
described in Bray et al. (1987). 

The Bray-Moss-Libby Model for the Aerothermochemistry 

The notion of bimodality extends to the the description of the velocity components in premixed turbulent 
combustion and leads to the Bray-Moss-Libby (BML) model of the aerothermochemistry of such 
combustion. To appreciate this extension consider the bivariate pdf 

P(u it c;x ) = a(x) 5(c) P(«, , 0,x ) + (J(*) 5(1 - c) P(m, , lyc) + 

y(x)f(u h c\x) (7) 

Figure 2 shows schematically such a joint pdf. If / (u, , c \x) is normalized, Eq. (4) again applies and with 

the approximation y <sc 1 the statistical behavior of u, and c are again dominated by contributions from 
reactants and products. Note that two of the three pdfs on the right side of Eq. (7) are conditional and 
describe the velocity within reactants and products. It is possible with current diagnostic techniques to 
measure with good spatial resolution the temperature and one or more velocity components and thus to 
determine experimentally the conditional pdfs appearing in Eq. (7) (cf. Moss 1980 and Shepherd and 
Cheng 1988). 

Equation (7) leads to significant results if y« 1. To illustrate consider the mean unconditional 
velocity component u t (x) which is given by Eqs. (5) and (7) as 

a i (x) = (l-d)a ir +d Uy (8) 

where u ir and u tp are the conditional mean values of the i-th velocity component within reactants and 
products respectively. We show these values in Fig. 2. Precise definitions of these quantities are as fol- 
lows: 

OO 

u ir (x)= l dui Ui P (Ui , 0;x) 
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«,>(*) = J du i u i P (“i > !pf) 

— oo 


( 9 ) 


According to the BM model the mean turbulent flux of all state variables in the i-th coordinate 
direction can be determined from the corresponding flux of the progress variable but from Eq. (7) we 
have 


pu/'c" 

— = — = C (1 - c) (u ip - u ir ) 
p 

Now a gradient model for this mean flux yields 


P <c” dc 

= - -v T — 

p OXi 


( 10 ) 


( 11 ) 


where v r is a positive turbulent exchange coefficient. If in connection with Fig. 1 we let x, correspond to 
x , then dcldx > 0 implying that pu ”c " < 0 and if gradient transport applies, then u p <u r . But we know 
from our earlier discussion that the velocity in the product stream downstream of the flame where c = 1 is 
(1 + x) times greater than that in the reactant stream. Given this overall increase in mean velocities it 
does not seem reasonable to expect the conditional velocities in products within the flame to be less than 
those within reactants. Accordingly, we see a potential difficulty with a gradient model for premixed tur- 
bulent flames if, as is the case in flames of practical interest, x » 1. 

In Bray et al. (1981), Masuya and Libby (1981) and Libby (1985) the BML model with all gradient 
assumptions avoided is applied to normal and oblique turbulent flames in premixed systems in order to 
clarify the applicability of the gradient model therein. In the next sections we sketch the analyses 
involved and review the results from them 
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APPLICATION TO NORMAL PREMIXED TURBULENT FLAMES 


We now consider application of the notions advanced in the previous section to the flame of Fig. 1. The 
equations needed for the analysis are as follows: 


-j- (P « 2 + P« 
dx 


" 2 \ _ _ dfL 

’ dx 


— (p u c + pu "c ") = w 
dx 


( 12 ) 


and 


d_ 

dx 

d_ 

dx 




p U — + pU " 3 


+ 2p«" 2 -^ = -2«"^--X u 


dx 


dx 


(13) 


, p« c 
U + 


pu" 2 c" 


du 




V = _p7dlL + -, 


+ p«"c" +p U"“ — =- c' -f- + « w - Xuc 

K dx dx dx 


The usual notation applies to Eqs. (12) and (13) which represent the first and second moment equations 
respectively for one dimensional flames. The quantities x« tuid describe dissipation effects which for 
the flames under consideration are due solely to chemical reaction, i.e., are proportional to w . Similarly 
the velocity-chemical reaction term u"w can be convincingly described by application of the BML 
model; it is found to be proportional to w, models for which are presently under development (cf. Bray et 
al. 1987) but which are not needed for present purposes* 

Some analysis permits Eqs. (12) and (13) to be reduced to two equations with c as the independent 
variable and with a dimensionless velocity intensity / = pu"^/p r u a and a dimensionless turbulent flux 
P = pu" c "/p r u Q as the two dependent variables. In the present discussion the latter variable is of more 
significance. In this formulation w is needed only if subsequent to finding the solutions for 1(c) and 
F(c) the spatial distributions are sought. However, comparison between theory and experiment is 

* II is worth noting that in these equations the effects of pressure fluctuations are neglected. In constant density turbulence models for such ef- 
fects have been painfully and carefully developed over a period of many years but these models do not apply to reacting flows with heat release 
since combustion induced pressure fluctuations result from volumetric sources, i.e., from an entirely different mechanism. Appropriate models 
for such fluctuations are not presently available. 
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conveniently carried out in terms of c so that generally the absence of a model for vv is not a shortcoming 
relative to the analysis of one dimensional flames. For premixed turbulent flows in two and three dimen- 
sions such a model is needed and research to that end is underway (cf. Bray etal. 1987). 

In developing the final equations for /(c) and F(c) the first of Eqs. (2) and (12) permit u to be 
eliminated while the second and third of Eqs. (12) are used to eliminate dp/d ; c and vv respectively. Clo- 
sure requires the third moment quantities pu" 3 and p u" 2 c", u", c" and the dissipation terms to be 
expressed in terms of the two dependent and the independent variables but with the exception of some 
inessential uncertainties in the conditional velocity statistics suitable models for these quantities free of 
gradient assumptions can be developed. For example, we have 




c (1 -g) 

1 + 1 C 


(14) 


In the present discussion two terms in Eqs. (13) are of particular interest; we refer to those involving 
dp/dx. If Eqs. (12) and (13) are specialized to constant density flows, e.g., by setting 1 = 0, then from 
Eqs. (14) we see that these terms vanish and indeed the reduction to standard equations for such flows is 
complete. The implication to be drawn from this consideration is that these pressure gradient terms 
account for effects operative in variable density turbulence, in particular an interaction between a force 
field associated with the pressure gradient and density fluctuations associated with heat release. In 
premixed turbulent flames this interaction is present only within the flame proper since both the force 
field and the density fluctuations are absent in the reactant and product streams on each side of the flame. 

Additional comments on this interaction are indicated. If conventional rather than Favre averaging 
is used, the interaction is contained within the resulting equations but is obscured by the clutter of a large 
number of terms involving density fluctuations. In more general variable density flows models for the 
multipliers of the several pressure gradient terms must be introduced as indicated for nonpremixed 
combustion by Jones (1980) and for high speed flows by Wilcox and Rubesin (1980). 

The final equations for 1(c) and F(c) involve singularities at c =0,1 as is to be expected since 
these points correspond to points at infinity in the x-variable but appropriate series expansions permit the 
numerical solutions to be initiated in the neighborhood of these end points. A shortcoming of the theory 



is that l a = u ' 2 0 lUg where u' 2 a is the intensity in the reactant stream upstream of the flame must be 
imposed, i.e., the turbulent flame speed is not calculated. This shortcoming can be turned into a virtue for 
the purpose of studying the structure of premixed turbulent flames since it can be argued that selection of 
1 0 so as to achieve agreement with experimental results assures that the predicted flame structure 
corresponds to realistic flames. In our original work we set I Q = 0.22 unless we had specific reasons to do 
otherwise but in recent years our confidence in this value has been reduced by evidence that a wide 
variety of values is found depending on the geometry of the flames and on the flameholding mechanism 
employed. The detailed reasons for this ambiguity are unknown.* 

Numerical Results for Normal Flames 

In Fig. 3 we show the distribution of the turbulent flux in terms of F(c) from both theory and from the 
experiment of Moss (1980). In the open flame studied by Moss x = 6.5 and 1„ =0.16. The two curves 
represent theoretical predictions based on slightly different models for the conditional statistics, differ- 
ences which are irrelevant for the present discussion. If it is recalled that gradient transport for these 
flames implies F < 0, we see that the entire flame structure exhibits countergradient diffusion. The expla- 
nation for such diffusion resides in the interaction between the pressure drop across the flame and the den- 
sity fluctuations which from Eq. (6) are seen to involve a relative intensity greater than unity, roughly 1.4. 
The pressure drop tends to accelerate the light parcels of product relative to the heavy parcels of reactants 
in the direction of high product concentration, a differential effect contrary to gradient transport. 

Several further comments are indicated. Calculations with the pressure gradient terms in Eqs. (13) 
suppressed yield F < 0 for all values of x. Moreover, as x — » 0 the theory predicts gradient transport and 
thus the expected behavior for nearly constant density turbulent flows. These results lend credibility to 
the explanation of interaction as the cause of countergradient diffusion and to the validity of the theory in 
general. The notion of countergradient transport in premixed turbulent flames is now widely accepted 


* A recent study (Libby 1987) removes the shortcoming with respect to 1 0 by invoking the Hakberg-Gosman condition and discusses in some de- 
tail the uncertainties in its value. 
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and has been observed in a variety of laboratory flames with simple geometries. Later we discuss its 
more general applicability. 
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THE OBLIQUE PREMIXED TURBULENT FLAME WITH CONSTRAINED MEAN STREAM- 
LINES 


In Fig. 4 we show schematically a premixed turbulent flame which is oblique to the reactant stream. Such 
flames can be established by a flameholder, e.g., a cylinder, in a duct carrying reactants. The classic 
experiment of this nature is due to Wright and Zukoski (1962) while the theoretical description of these 
flames based on the BML model is given by Masuya and Libby (1981). If we assume that the duct prohi- 
bits significant deflection of the mean streamlines both upstream and downstream of the flame, it is rea- 
sonable to consider as an idealization an infinite planar flame held at a specified angle 0 with undeflected 
mean streamlines. In this case description of the flow involves an analysis identical with that for purely 
normal flames and a second, subsequent analysis of the tangential velocity component which involves 
explicitly the flame angle 0. The treatment of the tangential flow requires determination of the intensity 
of the fluctuations of the v-velocity component and the tangential flux of the progress variable which in 
dimensionless form are /„ =pv" 2 /p r « 0 2 and F v = pv"c'7p r u a respectively. This latter variable is of 
principal concern in the context of the present discussion and from Eq.) is given by 

= c (1 -c) (v p -v r ) (15) 

P 


where v r and v p are the conditional tangential velocity components within reactants and products respec- 
tively. 


Gradient transport indicates that this mean flux is zero since the tangential gradient of the progress 
variable is zero and thus that the two conditional velocities in the tangential direction are equal. The 
implication from this result is that parcels of reactants and products have streamlines which differ only by 
the differences in the normal conditional velocities. However, there is a tangential force field in these 


flames arising from the x-wise gradient of the Reynolds shear stress p«"v"; the existence of this forces 
field can be seen from the mean conservation equation for tangential momentum which yields 


d — ttt, P r u o dc 

— p U V = — T — — 

dx tan 0 dx 


( 16 ) 


Note that when the heat release and thus the density fluctuations vanish, this force field is absent. 


This discussion establishes that oblique turbulent flames in premixed systems involve a tangential 
force field and density fluctuations. Thus according to our previous argument we can expect nongradient 
turbulent transport in the tangential direction. Indeed calculations show that v p >v r , that pv"c" > 0 and 
that the parcels of reactants and products follow different mean streamlines with the former exhibiting 
only small tangential velocity while the latter possesses large tangential velocities. This behavior is 
shown schematically in Fig. 4. We thus have an example of nongradient transport to add to the previously 
discussed case of countergradient transport. To date there have been no detailed measurements in oblique 
flames to assess the validity of this theory. 
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CONCLUDING REMARKS 


We show that when mean force fields and density fluctuations coexist, an interaction leads to turbulent 
transport which is not described by the usual gradient model. In the simple flow configurations associated 
with normal and oblique planar flames, the latter with constrained mean streamlines, countergradient and 
nongradient transport exists. There is no question that the notions suggested by these findings are con- 
ceptually important. Moreover, the experimental results of Heitor et al. (1987) relative to the complex 
flow associated with a baffle stabilized premixed turbulent flame establish that nongradient transport 
exists in complex flow configurations, i.e., in flows of applied interest. Within the context of the present 
discussion the first few sentences of this paper are worth quoting: 

In turbulent, premixed flames there arise source terms, explicitly set out below, in the conservation 
equations for the turbulent heat transfer rate and stresses that have no counterparts in non-reacting 
flows. Analysis (Masuya and Libby 1981; Bray, Libby and Moss 1985; Libby 1985) has shown that 
at least in the two idealized extremes with the flame either normal or oblique to the approaching 
reactants, and at practically important levels of heat release, these terms are sufficiently large to 
cause non-gradient transport of turbulent heat flux. This finding is important because it casts doubt 
on the applicability of turbulence models that use gradient-transport hypothesis. 

From a detailed study of the velocity and temperature fields Heitor et al. conclude that the interaction 
between pressure gradients and density fluctuations results in the largest contribution to the balance of the 
turbulent heat flux and that that flux is not aligned with the mean temperature gradient. 

The importance of nongradient transport on general variable density turbulence remains to be esta- 
blished. With respect to premixed turbulent flames we note that the question arises as to the influence of 
these processes arising within their structure on their orientation between reactant and product streams 
remains to be clarified. Answering this question requires further theoretical and computational efforts.* 

* If the provisional model for W ( X ) set forth in Bray et at. (1987) is validated, the effort required is largely computational. 
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Additional research is also required to determine the importance of the interaction under discussion 
in high speed flows. It would seem obvious that in scramjets with their large gradients of mean pressure 
and mean turbulent shear stresses this interaction must be operative but whether it plays an important role 
in the mixing and chemical processes in the flow is unclear. Certainly the present casual, uncritical exten- 
sion to supersonic chemically reacting flows of the phenomenology of constant density turbulence should 
raise healthy skepticism concerning the validity of the resulting predictions but unfortunately that skepti- 
cism does not appear to be shared by authors, reviewers and editors of current journal articles.** Similar 
uncertainties would appear to prevail for turbulent boundary layers in high speed flows so that additional 
research is indicated. 


** Occasionally there is a mild rejoinder warning the reader of potential fundamental difficulties. An example from a recent report is as follows: 
"The turbulence model we we have used is the k — £ model as described in ... Although Reference ... specifically addresses the question of 
compressible flows, it should be noted that turbulence models for compressible flows are not well-developed. Hence, the choice of the standard 
k — £ model in this situation cannot be regarded as definite." 
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